The genomic and evolutionary landscapes of anaplastic thyroid carcinoma.

Anaplastic thyroid carcinoma is arguably the most lethal human malignancy. It often co-occurs with differentiated thyroid cancers, yet the molecular origins of its aggressivity are unknown. We sequenced tumor DNA from 329 regions of thyroid cancer, including 213 from patients with primary anaplastic thyroid carcinomas. We also whole genome sequenced 9 patients using multi-region sequencing of both differentiated and anaplastic thyroid cancer components. Using these data, we demonstrate thatanaplastic thyroid carcinomas have a higher burden of mutations than other thyroid cancers, with distinct mutational signatures and molecular subtypes. Further, different cancer driver genes are mutated in anaplastic and differentiated thyroid carcinomas, even those arising in a single patient. Finally, we unambiguously demonstrate that anaplastic thyroid carcinomas share a genomic origin with co-occurring differentiated carcinomas and emerge from a common malignant field through acquisition of characteristic clonal driver mutations.

Early Cancer Detection in Li-Fraumeni Syndrome with Cell-Free DNA.

People with Li-Fraumeni syndrome (LFS) harbor a germline pathogenic variant in the TP53 tumor suppressor gene, face a near 100% lifetime risk of cancer, and routinely undergo intensive surveillance protocols. Liquid biopsy has become an attractive tool for a range of clinical applications, including early cancer detection. Here, we provide a proof-of-principle for a multimodal liquid biopsy assay that integrates a targeted gene panel, shallow whole-genome, and cell-free methylated DNA immunoprecipitation sequencing for the early detection of cancer in a longitudinal cohort of 89 LFS patients. Multimodal analysis increased our detection rate in patients with an active cancer diagnosis over uni-modal analysis and was able to detect cancer-associated signal(s) in carriers prior to diagnosis with conventional screening (positive predictive value = 67.6%, negative predictive value = 96.5%). Although adoption of liquid biopsy into current surveillance will require further clinical validation, this study provides a framework for individuals with LFS. SIGNIFICANCE: By utilizing an integrated cell-free DNA approach, liquid biopsy shows earlier detection of cancer in patients with LFS compared with current clinical surveillance methods such as imaging. Liquid biopsy provides improved accessibility and sensitivity, complementing current clinical surveillance methods to provide better care for these patients. See related commentary by Latham et al., p. 23. This article is featured in Selected Articles from This Issue, p. 5.

Loss of LRP1B expression drives acquired chemo and radio-resistance in HPV-positive head and neck cancer.

OBJECTIVES: Although human papillomavirus positive (HPV+) oropharyngeal squamous cell carcinoma (OPSCC) patients typically experience excellent survival, 15-20 % of patients recur after treatment with chemotherapy and radiation. Therefore, there is a need for biomarkers of treatment failure to guide treatment intensity. MATERIALS AND METHODS: Whole genome sequencing was carried out on HPV+OPSCC patients who were primarily treated with concurrent chemotherapy (cisplatin) and radiation. We then explored whether the loss of LRP1Bwas sufficient to drive an aggressive phenotype, and promote a resistance to cisplatin and radiation therapy both in vitro using HPV+ cell lines (93VU147T, UMSCC47, UWO37 and UWO23) and in vivo. RESULTS: Through integrative genomic analysis of three HPV+OPSCC tumour datasets, we identified that deletion of LRP1B was enriched in samples that recurred following chemo-radiation. Knockdown using siRNA in four HPV+ cell lines (UWO23, UWO37, UMSCC47 and 93VU147T) resulted in increased proliferation of all cases. CRISPR/Cas9 deletion of LRP1B in the same cell line panel demonstrated increased proliferation, clonogenic growth and migration, as well as resistance to both cisplatin and radiation in LRP1B deleted cells compared to their respective non-targeting control cells. Cell line derived xenograft studies indicated that the LRP1B knockout tumours were more resistant to cisplatin and radiation therapy compared to their controls invivo. CONCLUSION: Taken together, our work implicates LRP1B deletion as a potential biomarker for identifying treatment resistant HPV+ OPSCC cases.

Identifying Mechanisms of Resistance by Circulating Tumor DNA in EVOLVE, a Phase II Trial of Cediranib Plus Olaparib for Ovarian Cancer at Time of PARP Inhibitor Progression.

PURPOSE: To evaluate the use of blood cell-free DNA (cfDNA) to identify emerging mechanisms of resistance to PARP inhibitors (PARPi) in high-grade serous ovarian cancer (HGSOC). EXPERIMENTAL DESIGN: We used targeted sequencing (TS) to analyze 78 longitudinal cfDNA samples collected from 30 patients with HGSOC enrolled in a phase II clinical trial evaluating cediranib (VEGF inhibitor) plus olaparib (PARPi) after progression on PARPi alone. cfDNA was collected at baseline, before treatment cycle 2, and at end of treatment. These were compared with whole-exome sequencing (WES) of baseline tumor tissues. RESULTS: At baseline (time of initial PARPi progression), cfDNA tumor fractions were 0.2% to 67% (median, 3.25%), and patients with high ctDNA levels (>15%) had a higher tumor burden (sum of target lesions; P = 0.043). Across all timepoints, cfDNA detected 74.4% of mutations known from prior tumor WES, including three of five expected BRCA1/2 reversion mutations. In addition, cfDNA identified 10 novel mutations not detected by WES, including seven TP53 mutations annotated as pathogenic by ClinVar. cfDNA fragmentation analysis attributed five of these novel TP53 mutations to clonal hematopoiesis of indeterminate potential (CHIP). At baseline, samples with significant differences in mutant fragment size distribution had shorter time to progression (P = 0.001). CONCLUSIONS: Longitudinal testing of cfDNA by TS provides a noninvasive tool for detection of tumor-derived mutations and mechanisms of PARPi resistance that may aid in directing patients to appropriate therapeutic strategies. With cfDNA fragmentation analyses, CHIP was identified in several patients and warrants further investigation.

Serological abnormalities that predict progression to systemic autoimmune rheumatic diseases in antinuclear antibody-positive individuals.

OBJECTIVE: We investigated the autoantibody (autoAb) profiles in ANA+ individuals lacking systemic autoimmune rheumatic disease (SARD) and early SARD patients to determine the key differences between these groups and identify factors that are associated with an increased risk of symptomatic progression within the next 2 years in ANA+ individuals. METHODS: Using custom antigen (Ag) microarrays, 144 IgM and IgG autoAbs were surveyed in 84 asymptomatic and 123 symptomatic (48 UCTD and 75 SARD patients) ANA+ individuals. AutoAbs were compared in ANA+ individuals lacking a SARD diagnosis with ≥2 years follow-up (n = 52), including all those who demonstrated progression (n = 14) during this period, with changes over time assessed in a representative subset. RESULTS: We show that ANA+ individuals have autoAb to many self-Ags that are not being captured by current screening techniques and very high levels of these autoAbs are predominantly restricted to early SARD patients, with SLE patients displaying reactivity to many more autoAgs than the other groups. In general, the symptoms that developed in progressors mirrored those seen in SARD patients with similar patterns of autoAbs. Only anti-Ro52 Abs were found to predict progression (positive predictive value 46%, negative predictive value 89%). Surprisingly, over 2 years of follow-up the levels of autoAbs remained remarkably stable regardless of whether individuals progressed or not. CONCLUSION: Our findings strongly argue that development of assays with an expanded set of auto-Ags and enhanced dynamic range would improve the diagnostic and prognostic ability of autoAb testing.

Flavopiridol causes cell cycle inhibition and demonstrates anti-cancer activity in anaplastic thyroid cancer models.

Anaplastic thyroid cancer (ATC) is a rare, but nearly uniformly fatal disease that is typically resistant to chemotherapy and radiation. Alternative strategies to target this cancer at a molecular level are necessary in order to improve dismal outcomes for ATC patients. We examined the effects of flavopiridol, a CDK inhibitor, in a panel of ATC cell lines. When cell lines were treated over a ten-point concentration range, CAL62, KMH2 and BHT-101 cell lines had a sub micromolar half-maximal inhibitory concentration, while no effect was seen in the non-cancerous cell line IMR-90. Flavopiridol treatment resulted in decreased levels of the cell cycle proteins CDK9 and MCL1, and induced cell cycle arrest. Flavopiridol also decreased the in vitro ability of ATC cells to form colonies and impeded migration using a transwell migration assay. In vivo, flavopiridol decreased tumor weight and tumor volume over time in a patient-derived xenograft model of ATC. Given the observed in vitro and in vivo activity, flavopiridol warrants further investigation for treatment of ATC.

Sex differences in oncogenic mutational processes.

Sex differences have been observed in multiple facets of cancer epidemiology, treatment and biology, and in most cancers outside the sex organs. Efforts to link these clinical differences to specific molecular features have focused on somatic mutations within the coding regions of the genome. Here we report a pan-cancer analysis of sex differences in whole genomes of 1983 tumours of 28 subtypes as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium. We both confirm the results of exome studies, and also uncover previously undescribed sex differences. These include sex-biases in coding and non-coding cancer drivers, mutation prevalence and strikingly, in mutational signatures related to underlying mutational processes. These results underline the pervasiveness of molecular sex differences and strengthen the call for increased consideration of sex in molecular cancer research.

Transgenerational epigenetic and transcriptomic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin exposure in rat.

In rats, direct exposure to TCDD causes myriad toxicities. Exposed rats experience hepatotoxicity, wasting syndrome and immune suppression, amongst others. "Inherited exposure", as occurs in the F3 generation of directly exposed F0 animals, has also been shown to cause toxicity: both male and female F3 rats demonstrate an increased incidence of adult onset disease, females also display reproductive abnormalities and increased incidence of ovarian diseases while males show increased incidence of kidney disease and an altered sperm epigenome. Here, we explore the hepatic transcriptomic profile of male and female F3 Sprague-Dawley rats bred through the paternal germ line from F0 dams exposed to a single dose of TCDD (0, 30, 100, 300 or 1000 ng/kg body weight) by oral gavage. We hypothesize that RNA transcripts with altered abundance in livers of unexposed F3 progeny of treated F0 Sprague-Dawley rats may result from epigenetic modifications to the genome. We further survey patterns of differential methylation within male F3 rat testis. Female F3 rats demonstrated more TCDD-mediated hepatic transcriptomic changes than males, with differences primarily in the lowest dose group. In testis from male F3 rats, multiple olfactory receptors displayed patterns of differential methylation. Hypermethylation of Egfr and Mc5r among testes from TCDD lineage rats was observed, but without corresponding changes in hepatic mRNA abundance. Further studies examining these differences in other tissue types are warranted.

Sex disparities in head & neck cancer driver genes: An analysis of the TCGA dataset.

OBJECTIVES: Survival in head and neck squamous cell carcinoma (HNSCC) has been associated with patient sex, typically with males experiencing poorer outcomes. It is unclear if this disparity is based in divergent tumor biology. We analyzed the TCGA HNSCC cohort to uncover disparities in the somatic single nucleotide variation (SNV), copy number alteration (CNA) and mRNA abundance profiles between males and females. Critically, we stratified our results by tumor HPV status to control for this significant confounder. METHODS: SNV, CNA and mRNA abundance differences between males and females were compared separately for the HPV-positive (n = 67) and negative (n = 431) TCGA HNSCC cohorts. Overall survival outcomes were compared in males and females in both HPV-positive and HPV-negative subsets of patients. RESULTS: Females were found to have poorer overall survival than males (p = 0.048), largely due to higher rates of HPV-positive disease among men. SNV analysis revealed that in HPV-positive disease, there were no differences by sex after accounting for the false discovery rate (FDR). In HPV-negative tumors, BRWD3 mutations occurred more frequently in the tumors of female patients compared to males after adjusting for the FDR (p = 0.02). Further, HPV-negative BRWD3 mutant tumors were found to have significantly worse 5-year overall survival compared to wildtype on multivariate analysis (p = 0.02). There were 88 heterozygous deletions and 14 amplifications that were differentially altered between male and female HPV-negative tumors and associated with expression changes. Pathway analysis of these genes revealed that tumors from males were enriched in five pathways including chemokine and phosphophatidylinositol signaling. CONCLUSIONS: Reanalysis of the TCGA HNSCC dataset stratified by sex revealed that males in this cohort had a significant survival advantage, due to a higher proportion of HPV-positive disease. Mutations in BRWD3 were more frequent in HPV-negative tumors of females and were associated with poorer overall survival. BRWD3 may represent a novel biomarker of patient outcomes, but will require additional validation.

Genomic and human papillomavirus profiling of an oral cancer cohort identifies TP53 as a predictor of overall survival.

BACKGROUND: The genomic landscape of head and neck cancer has been reported through The Cancer Genome Atlas project. We attempt to determine if high-risk human papillomavirus (HPV) or frequently mutated genes are correlated with survival in an oral cancer cohort. METHODS: Patient demographic data along with data from final pathology was collected. Tumor DNA was analyzed using a custom Illumina targeted sequencing panel. Five high-risk HPV types were tested by qPCR. Statistical analyses were used to identify associations between patient outcome and mutational status. RESULTS: High-risk HPV types were identified in 7% of cases; HPV status was not associated with survival. Mutations were identified in TP53, TERT promoter, & PIK3CA. Mutations in TP53 were significantly associated with poorer overall survival on multi-variate analysis (p = 0.03). CONCLUSIONS: Mutations in TP53 were associated with poor patient survival. Expanding our sample size may identify further predictors of outcome to direct customized cancer care.

Rethinking Lupus Nephritis Classification on a Molecular Level.

The International Society of Nephrology/Renal Pathology Society (ISN/RPS) lupus nephritis (LN) classification is under reconsideration, given challenges with inter-rater reliability and resultant inconsistent relationship with treatment response. Integration of molecular classifiers into histologic evaluation can improve diagnostic precision and identify therapeutic targets. This study described the relationship between histological and molecular phenotypes and clinical responses in LN. Renal compartmental mRNA abundance was measured in 54 biopsy specimens from LN patients and correlated to ISN/RPS classification and individual histologic lesions. A subset of transcripts was also evaluated in sequential biopsies of a separate longitudinal cohort of 36 patients with paired samples obtained at the time of flare and at follow up. Unsupervised clustering based on mRNA abundance did not demonstrate a relationship with the (ISN/RPS) classification, nor did univariate statistical analysis. Exploratory analyses suggested a correlation with individual histologic lesions. Glomerular FN1 (fibronectin), SPP1 (secreted phosphoprotein 1), and LGALS3 (galectin 3) abundance correlated with disease activity and changed following treatment. Exploratory analyses suggested relationships between specific transcripts and individual histologic lesions, with the important representation of interferon-regulated genes. Our findings suggested that the current LN classification could be refined by the inclusion of molecular descriptors. Combining molecular and pathologic kidney biopsy phenotypes may hold promise to better classify disease and identify actionable treatment targets and merits further exploration in larger cohorts.

Comparative toxicoproteogenomics of mouse and rat liver identifies TCDD-resistance genes.

The aryl hydrocarbon receptor (AHR) mediates many toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). However, the AHR alone does not explain the widely different outcomes among organisms. To identify the other factors involved, we evaluated three transgenic mouse lines, each expressing a different rat AHR isoform (rWT, DEL, and INS) providing widely different resistance to TCDD toxicity, as well as C57BL/6 and DBA/2 mice which exhibit a ~ tenfold divergence in TCDD sensitivity (exposures of 5-1000 µg/kg TCDD). We supplement these with whole-genome sequencing, together with transcriptomic and proteomic analyses of the corresponding rat models, Long-Evans (L-E) and Han/Wistar (H/W) rats (having a ~ 1000-fold difference in their TCDD sensitivities; 100 µg/kg TCDD), to identify genes associated with TCDD-response phenotypes. Overall, we identified up to 50% of genes with altered mRNA abundance following TCDD exposure are associated with a single AHR isoform (33.8%, 11.7%, 5.2% and 0.3% of 3076 genes altered unique to rWT, DEL, C57BL/6 and INS respectively following 1000 µg/kg TCDD). Hepatic Pxdc1 was significantly repressed in all three TCDD-sensitive animal models (C57BL/6 and rWT mice, and L-E rat) after TCDD exposure. Three genes, including Cxxc5, Sugp1 and Hgfac, demonstrated different AHRE-1 (full) motif occurrences within their promoter regions between rat strains, as well as different patterns of mRNA abundance. Several hepatic proteins showed parallel up- or downward alterations with their RNAs, with three genes (SNRK, IGTP and IMPA2) showing consistent, strain-dependent changes. These data show the value of integrating genomic, transcriptomic and proteomic evidence across multi-species models in toxicologic studies.

2,3,7,8-Tetrachlorodibenzo-p-dioxin modifies alternative splicing in mouse liver.

Alternative splicing is a co-transcriptional mechanism that generates protein diversity by including or excluding exons in different combinations, thereby expanding the diversity of protein isoforms of a single gene. Abnormalities in this process can result in deleterious effects to human health, and several xenobiotics are known to interfere with splicing regulation through multiple mechanisms. These changes could lead to human diseases such as cancer, neurological disorders, autoimmune diseases, and developmental disorders. 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is an environmental contaminant generated as a byproduct of various industrial activities. Exposure to this dioxin has been linked to a wide range of pathologies through the alteration of multiple cellular processes. However, the effects of TCDD exposure on alternative splicing have not yet been studied. Here, we investigated whether a single po. dose of 5 µg/kg or 500 µg/kg TCDD influence hepatic alternative splicing in adult male C57BL/6Kou mouse. We identified several genes whose alternative splicing of precursor messenger RNAs was modified following TCDD exposure. In particular, we demonstrated that alternative splicing of Cyp1a1, Ahrr, and Actn1 was significantly altered after TCDD treatment. These findings show that the exposure to TCDD has an impact on alternative-splicing, and suggest a new avenue for understanding TCDD-mediated toxicity and pathogenesis.

A readout of metabolic efficiency in arylamine N-acetyltransferase-deficient mice reveals minor energy metabolism changes.

Recent studies have revealed a possible link between the activities of polymorphic arylamine N-acetyltransferases (NATs) and energy metabolism. We used a Nat1/Nat2 double knockout (KO) mouse model to demonstrate that ablation of the two Nat genes is associated with modest, intermittent alterations in respiratory exchange rate. Pyruvate tolerance tests show that double KO mice have attenuated hepatic gluconeogenesis when maintained on a high-fat/high-sucrose diet. Absence of the two Nat genes also leads to an increase in the hepatic concentration of coenzyme A in mice fed a high-fat/high-sucrose diet. Our results suggest a modest involvement of NAT in energy metabolism in mice, which is consistent with the absence of major phenotypic deregulation of energy metabolism in slow human acetylators.

Transcriptomic Impact of IMA-08401, a Novel AHR Agonist Resembling Laquinimod, on Rat Liver.

IMA-08401 (C2) is a novel aryl hydrocarbon receptor (AHR) agonist and selective AHR modulator (SAHRM) that is structurally similar to laquinimod (LAQ). Both compounds are converted to the AHR-active metabolite DELAQ (IMA-06201) in vivo. SAHRMs have been proposed as therapeutic options for various autoimmune disorders. Clinical trials on LAQ have not reported any significant toxic outcomes and C2 has shown low toxicity in rats; however, their functional resemblance to the highly toxic AHR agonist 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) raises questions. Here, we characterize the hepatic transcriptomic changes induced by acute (single-dose) and subacute exposure (repeated dosing for 5 days followed by a 5-day recovery period) to C2 in Sprague-Dawley rats. Exposure to C2 leads to activation of the AHR, as shown by altered transcription of Cyp1a1. We identify a heightened response early after exposure that drops off by day 10. Acute exposure to C2 leads to changes to transcription of genes involved in antiviral and antibacterial responses, which highlights the immunomodulator effects of this AHR agonist. Subacute exposure causes an oxidative stress response in the liver, the consequences of which require further study on target tissues such as the CNS and immune system, both of which may be compromised in this patient population.

Mutational analysis of head and neck squamous cell carcinoma stratified by smoking status.

Smoking has historically been recognized as a negative prognostic factor in head and neck squamous cell carcinoma (HNSCC). This study aimed to assess the mutational differences between heavy smokers (>20 pack years) and never smokers among the HNSCC patients within The Cancer Genome Atlas (TCGA). Single nucleotide variation and copy number aberration differences between heavy smokers and never smokers were compared within human papillomavirus-positive (HPV-positive) (n = 67) and HPV-negative (n = 431) TCGA cohorts with HNSCC, and the impact of these mutations on survival were assessed. No genes were differentially mutated between smoking and never-smoking patients with HPV-positive tumors. By contrast, in HPV-negative tumors, NSD1 and COL1A11 were found to be more frequently mutated in heavy smokers, while CASP8 was more frequently altered in never smokers. HPV-negative patients with NSD1 mutations experienced significantly improved overall survival compared with NSD1 WT patients. This improved prognosis was validated in an independent cohort of 77 oral cavity cancer patients and a meta-analysis that included 2 additional data sets (688 total patients, hazard ratio for death 0.44, 95% CI, 0.30-0.65). NSD1 mutations are more common in HPV-negative heavy smokers, define a cohort with favorable prognosis, and may represent a clinically useful biomarker to guide treatment deintensification for HPV-negative patients.

BPG: Seamless, automated and interactive visualization of scientific data.

BACKGROUND: We introduce BPG, a framework for generating publication-quality, highly-customizable plots in the R statistical environment. RESULTS: This open-source package includes multiple methods of displaying high-dimensional datasets and facilitates generation of complex multi-panel figures, making it suitable for complex datasets. A web-based interactive tool allows online figure customization, from which R code can be downloaded for integration with computational pipelines. CONCLUSION: BPG provides a new approach for linking interactive and scripted data visualization and is available at http://labs.oicr.on.ca/boutros-lab/software/bpg or via CRAN at https://cran.r-project.org/web/packages/BoutrosLab.plotting.general.

Lestaurtinib is a potent inhibitor of anaplastic thyroid cancer cell line models.

Anaplastic thyroid cancer (ATC) is a rare and lethal human malignancy with no known effective therapies in the majority of cases. Despite the use of conventional treatments such as chemotherapy, radiation and surgical resection, this disease remains almost universally fatal. In the present study, we identified the JAK2 inhibitor Lestaurtinib as a potent compound when testing against 13 ATC cell lines. Lestaurtinib demonstrated a potent antiproliferative effect in vitro at nanomolar concentrations. Furthermore, Lestaurtinib impeded cell migration and the ability to form colonies from single cells using scratch-wound and colony formation assays, respectively. Flow cytometry was used for cell cycle analysis following drug treatment and demonstrated arrest at the G2/M phase of the cell cycle, indicative of a cytostatic effect. In vivo studies using the chick chorioallantoic membrane xenograft models demonstrated that treatment with Lestaurtinib resulted in a significant decrease in endpoint tumor volume and vascularity using power Doppler ultrasound imaging. Overall, this study provides evidence that Lestaurtinib is a potent antiproliferative agent with potential antiangiogenic activity that warrants further investigation as a targeted therapy for ATC.

Identification of a neutrophil-related gene expression signature that is enriched in adult systemic lupus erythematosus patients with active nephritis: Clinical/pathologic associations and etiologic mechanisms.

Both a lack of biomarkers and relatively ineffective treatments constitute impediments to management of lupus nephritis (LN). Here we used gene expression microarrays to contrast the transcriptomic profiles of active SLE patients with and without LN to identify potential biomarkers for this condition. RNA isolated from whole peripheral blood of active SLE patients was used for transcriptomic profiling and the data analyzed by linear modeling, with corrections for multiple testing. Results were validated in a second cohort of SLE patients, using NanoString technology. The majority of genes demonstrating altered transcript abundance between patients with and without LN were neutrophil-related. Findings in the validation cohort confirmed this observation and showed that levels of RNA abundance in renal remission were similar to active patients without LN. In secondary analyses, RNA abundance correlated with disease activity, hematuria and proteinuria, but not renal biopsy changes. As abundance levels of the individual transcripts correlated strongly with each other, a composite neutrophil score was generated by summing all levels before examining additional correlations. There was a modest correlation between the neutrophil score and the blood neutrophil count, which was largely driven by the dose of glucocorticosteroids and not the proportion of low density and/or activated neutrophils. Analysis of longitudinal data revealed no correlation between baseline neutrophil score or changes over the first year of follow-up with subsequent renal flare or treatment outcomes, respectively. The findings argue that although the neutrophil score is associated with LN, its clinical utility as a biomarker may be limited.

Genomic hallmarks of localized, non-indolent prostate cancer.

Prostate tumours are highly variable in their response to therapies, but clinically available prognostic factors can explain only a fraction of this heterogeneity. Here we analysed 200 whole-genome sequences and 277 additional whole-exome sequences from localized, non-indolent prostate tumours with similar clinical risk profiles, and carried out RNA and methylation analyses in a subset. These tumours had a paucity of clinically actionable single nucleotide variants, unlike those seen in metastatic disease. Rather, a significant proportion of tumours harboured recurrent non-coding aberrations, large-scale genomic rearrangements, and alterations in which an inversion repressed transcription within its boundaries. Local hypermutation events were frequent, and correlated with specific genomic profiles. Numerous molecular aberrations were prognostic for disease recurrence, including several DNA methylation events, and a signature comprised of these aberrations outperformed well-described prognostic biomarkers. We suggest that intensified treatment of genomically aggressive localized prostate cancer may improve cure rates.